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We  established  a  three-dimensional  finite  element  model  of  thermoelectric  module  based  on 
low-temperature  thermoelectric  material  bismuth  telluride  and  medium-temperature  thermoelectric 
material  filled-skutterudite.  The  material  properties  of  the  thermoelectric  materials  such  as  the  Seebeck 
coefficient,  thermal  conductivity,  and  electrical  conductivity  are  temperature  dependent.  Based  on  the 
formal  model,  multi-stage  models  consist  of  low-  and  medium-temperature  thermoelectric  modules 
are  proposed.  The  effect  of  input  energy  on  performance  of  solar  thermoelectric  generator  is  considered 
according  to  the  real  operating  condition.  Results  show  that,  reasonable  thermal  design  of  solar  thermo¬ 
electric  generator  can  take  full  advantage  of  the  characteristics  of  thermoelectric  materials  and  effectively 
improve  the  performance  of  power  generation. 
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1.  Introduction 

Solar  thermoelectric  generator,  which  uses  concentrated  solar 
radiation  as  heat  source,  has  been  studied  for  many  years.  Chen 
developed  a  thermodynamic  model  for  discussing  the  perfor¬ 
mance  of  solar-driven  thermoelectric  generator  [1].  Scherrer 
et  al.  developed  a  numerical  model  of  solar  thermoelectric 
generator  base  on  skutterudite  for  satellite  [2].  Omer  and  Infield 
presented  an  improved  theoretical  model  of  solar  thermoelectric 
generator  which  includes  all  parameters  that  have  an  influence 
on  the  heat  transfer  process,  and  the  modeling  results  are 
compared  with  experimental  data  of  the  commercial 
thermoelectric  module  in  power  generation  mode  [3].  Li  et  al. 
presented  a  prototype  concentration  solar  thermoelectric  gener¬ 
ator  (CTG)  and  a  discrete  numerical  model  for  the  evaluation  of 
the  whole  system  [4].  All  these  theoretical  and  experimental 
studies  lay  a  foundation  for  the  development  of  solar  thermo¬ 
electric  generator. 

There  are  many  advantages  of  solar  thermoelectric  generator, 
such  as  endless  shelf  life,  simple  structure,  no  moving  parts,  silent 
in  operation  and  no  pollution.  However,  the  low  conversion  effi¬ 
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ciency  makes  it  hard  to  be  widely  used.  The  conversion  efficiency 
of  solar  thermoelectric  generator  is  mainly  restricted  by  thermo¬ 
electric  materials.  Over  the  past  several  years,  a  number  of  high 
performance  thermoelectric  materials  have  been  developed  [5,6], 
and  some  of  them  are  available  commercially.  Bismuth  telluride 
is  a  favorable  low-temperature  thermoelectric  material,  which  ex¬ 
erts  better  property  at  low  temperature  (25-225  °C)  [7,8].  While 
filled-skutterudite  is  a  good  medium-temperature  thermoelectric 
material,  that  can  be  used  for  wide  temperature  range  generation 
(25-525  °C)  [9]. 

In  addition  to  the  research  of  thermoelectric  materials,  reason¬ 
able  thermal  design  and  management  of  thermoelectric  generator 
is  equally  important  for  improving  the  generating  performance. 
Chen  et  al.  established  a  theoretical  model  of  a  two-stage  semicon¬ 
ductor  thermoelectric  generator  [10].  El-Genk  and  Saber  presented 
a  1-D  analytical  model  of  segmented  thermoelectric  unicouple  for 
operation  between  973  and  300  K  [11].  In  the  present  paper,  we 
construct  a  three-dimensional  finite  element  model  [12,13]  for  dis¬ 
cussing  the  performance  characteristics  of  low-  and  medium- 
temperature  thermoelectric  modules.  The  material  properties  of 
the  thermoelectric  materials  such  as  the  Seebeck  coefficient,  ther¬ 
mal  conductivity,  and  electrical  conductivity  are  temperature 
dependent.  Based  on  the  singe-stage  models,  we  present  two-  and 
three-stage  models  of  the  solar  thermoelectric  generator  to  discuss 
its  performance  characteristics.  The  effect  of  input  energy  on  perfor¬ 
mance  of  solar  thermoelectric  generator  is  considered  according  to 
the  real  operating  condition. 
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Nomenclature 

C 

specific  heat  capacity  (J/kg  K) 

Greek  symbols 

T 

temperature  (°C) 

p 

density  (kg/m3) 

Q 

heat  generation  rate  per  unit  volume  (W/m3) 

w 

thermal  conductivity  matrix  (W/mK) 

3 

heat  flux  vector  (W/m2) 

M 

electrical  conductivity  matrix,  S/m 

J 

electric  current  density  vector  (A/m2) 

[a] 

Seebeck  coefficient  matrix  (V/K) 

D 

electric  flux  density  vector  (C/m2) 

M 

dielectric  permittivity  matrix  (F/m) 

E 

electric  field  intensity  vector  (V/m) 

(P 

electric  scalar  potential  (V) 

I 

load  current  (A) 

hSTEG 

total  efficiency  of  solar  thermoelectric  generator 

Rl 

electric  load  (Q) 

Pte 

thermo  electric  conversion  efficiency 

P 

output  power  (W) 

hopt 

optical  focusing  efficiency 

Qin 

input  power  (W) 

ha 

absorptivity  of  the  heat  collector 

qsolar 

solar  radiation  flux  density,  1000  W/m2 

Q 

total  concentration  ratio  of  optical  focusing  system 

2.  Mathematical  model  of  solar  thermoelectric  generator  2.2.  Model  parameters  of  solar  thermoelectric  generator 


2A.  Governing  equations  of  coupled  thermoelectricity 

The  heat  flow  equation  in  thermoelectric  analysis  can  be  ex¬ 
pressed  as: 

pcg  +  V-q  =  q  (1) 

where  p ,  C,  T,  q  and  q  stand  for  density,  specific  heat  capacity,  tem¬ 
perature,  heat  generation  rate  per  unit  volume  and  heat  flux  vector, 
respectively. 

The  continuity  equations  of  electric  charge  is 

V.  (/+!). °  (2) 

where  J  is  the  electric  current  density  vector,  D  is  the  electric  flux 
density  vector. 

Above  two  equations  are  coupled  by  the  set  of  thermoelectric 
constitutive  equations  [12], 

q  =  T[a] -J-  [A]  -  VT  (3) 

J=[a]-(£-[a]-VT)  (4) 

where  [2]  is  the  thermal  conductivity  matrix,  [cr]  is  the  electrical 
conductivity  matrix,  and  [a]  is  the  Seebeck  coefficient  matrix. 
E  =  -Vcp  is  the  electric  field  intensity  vector,  where  cp  is  the  electric 
scalar  potential. 

The  coupled  equations  of  thermoelectricity  can  be  obtained 


from  above  equations, 

pcg  +  V.(T[a].J)-V.([A].VT)  =  q 

(5) 

V  '  ([£] '  V^)  +  V  ' M  '  VT)  +  V  ' '  V<?)  =  ° 

(6) 

where  [s]  is  the  dielectric  permittivity  matrix. 

In  the  present  steady-state  model,  material  properties  of  all 
components  are  considered  to  be  isotropic.  The  coupled  equations 
of  thermoelectricity  can  be  written  as  follows: 

V  •  ( TotJ )  -  V  •  (AVT)  =  q 

(7) 

V  •  ( aocVT )  +  V  ■  (<rV<p)  -  0 

(8) 

The  equations  described  in  the  previous  section  are  solved  using 
the  ANSYS  which  is  the  commercial  software  based  on  the  finite 
element  method  and  is  applicable  to  multiphysics  problems  [13]. 


2.2 A.  Geometric  model 

A  typical  solar  thermoelectric  generator  consists  of  an  optical 
focusing  system,  heat  collector,  thermoelectric  module,  and  cool¬ 
ing  plate.  The  thermoelectric  module  is  the  major  component  of 
the  solar  thermoelectric  generator.  As  shown  in  Fig.  la,  the 
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Fig.  1.  Geometric  model  of  thermoelectric  module  and  thermoelectric  unicouple. 
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Table  1 

Geometric  size  of  thermoelectric  unicouple. 


Unicouple 

Unicouple  section 

Conductive  film 

Unicouple 

length  (mm) 

(mm  x  mm) 

thickness  (mm) 

gap  (mm) 

1.8 

1  x  1 

0.1 

0.4 

thermoelectric  module  consists  of  conductive  film,  thermal  insu¬ 
lating  material  and  98  pairs  of  unicouples.  In  order  to  facilitate 
the  calculation,  we  analyzed  only  one  unicouple  in  the  present  pa¬ 
per.  The  geometric  model  of  the  unicouple  is  shown  in  Fig.  lb.  The 
p-  and  n-type  thermoelectric  materials  are  connected  by  the  con¬ 
ductive  film,  and  the  sawtooth  line  stand  for  electric  load,  the  geo¬ 
metric  size  of  the  unicouple  is  shown  in  Table  1. 

2.2.2.  Material  properties 

The  selection  of  thermoelectric  materials  directly  affects  the 
performance  of  the  generator.  We  used  bismuth  telluride  and 
filled-skutterudite  as  materials  of  the  low-  and  medium-tempera¬ 
ture  modules,  respectively.  As  the  solar  thermoelectric  generator  is 
operated  under  large  temperature  differences,  the  material  proper¬ 
ties  of  the  thermoelectric  materials  such  as  the  Seebeck  coefficient, 
thermal  conductivity,  and  electrical  conductivity  are  temperature 
dependent  [7-9].  The  thermal  conductivity  of  the  electrically  con¬ 
ductive  film  is  300  W/m  K.  The  contact  resistance  is  not  taken  into 
account  in  the  present  model. 

2.2.3.  Boundary  conditions 

The  solar  radiation  focused  by  optical  focusing  system  is  ab¬ 
sorbed  by  the  heat  collector  for  conversion  into  heat  energy.  The 
input  heat  flux  density  can  be  determined  as  follows  [4]: 

Qin  —  Qsolar^-gQoptQ a  (^) 

where  qso/ar  is  the  solar  radiation  flux  density  ( qsoiar=  1000  W/m2) 
[1,14],  Cg  is  the  total  concentration  ratio  of  the  optical  focusing 
system,  popt  is  the  optical  focusing  efficiency,  and  r\a  is  the  absorp¬ 
tivity  of  the  heat  collector. 

The  actual  generator  is  cooled  by  cooling  water.  In  the  present 
paper,  the  cooling  of  generator  is  assumed  to  be  very  good  and 
the  first  boundary  condition  is  applied  to  the  cold  end  of  the  ther¬ 
moelectric  module  with  a  fixed  value  of  25  °C. 

In  operation,  heat  losses  mainly  include  convection  loss  and 
radiation  loss.  Convection  heat  loss  refers  to  the  convection  be¬ 
tween  the  surface  of  the  generator  and  ambient  air.  Radiation  heat 
loss  is  the  radiation  of  the  surface  to  ambient.  Since  the  main  work 
of  this  paper  is  to  study  the  performance  of  the  solar  thermoelec¬ 
tric  generator  under  ideal  conditions,  the  heat  losses  due  to  con¬ 
vection  and  radiation  are  not  taken  into  account. 

2.3.  Output  power  and  conversion  efficiency  of  solar  thermoelectric 
generator 

The  output  power  is  defined  as  follows, 

P  =  I2Rl  (10) 

where  I  and  RL  stand  for  load  current  and  electric  load,  respectively. 
So  the  thermoelectric  conversion  efficiency  can  be  obtained  by, 

QtE  =  P/Qin  (11) 

where  Qin  stands  for  input  power.  The  total  efficiency  of  solar  ther¬ 
moelectric  generator  is  the  product  of  the  thermoelectric  conver¬ 
sion  efficiency  rjTE ,  the  optical  focusing  efficiency  rjopU  and  the 
absorptivity  of  the  heat  collector  r\a. 

QsTEG  =  VTEVoptVa  (12) 


3.  Results  and  discussion 

3.1.  Single-stage  thermoelectric  modules 

The  performance  of  the  single-stage  thermoelectric  modules  di¬ 
rectly  affects  the  overall  solar  utilization  of  solar  thermoelectric 
generator.  Therefore,  the  analysis  of  generation  performances  and 
characteristics  for  the  low-  and  medium-temperature  thermoelec¬ 
tric  modules  should  be  carried  out  prior  to  the  thermal  design  of 
the  solar  thermoelectric  generator. 

Figs.  2  and  3  show  that  the  conversion  efficiency  and  output 
power  of  each  module  rise  with  increasing  temperature  difference. 
When  the  temperature  difference  rises  from  0  °C  to  200  °C,  the 
conversion  efficiency  of  the  bismuth  telluride  unicouple  increases 
from  0%  to  10.82%  and  the  output  power  from  0W  to  0.0279  W, 
respectively.  Meanwhile,  the  conversion  efficiency  and  output 
power  of  the  skutterudite  unicouple  reaches  9.24%  and  0.1009  W, 
respectively,  at  the  temperature  difference  of  500  °C.  It  also  can 
be  seen  that  the  performance  of  the  bismuth  telluride  unicouple 
is  better  than  that  of  the  skutterudite  one  when  the  unicouples 
are  operating  between  25  °C  and  225  °C.  Therefore,  effective 
enhancement  of  conversion  efficiency  could  be  achieved  by  using 
multi-stage  thermoelectric  module  which  can  make  full  use  of 
the  advantages  of  each  thermoelectric  material. 

3.2.  Two-stage  thermoelectric  module 

In  above,  we  analyze  the  performance  of  the  low-  and  medium- 
temperature  thermoelectric  unicouples.  The  analysis  is  the  founda¬ 
tion  for  researching  the  two-stage  thermoelectric  module,  which 
consists  of  a  low-temperature  module  and  a  medium-temperature 
module.  The  medium-temperature  module  is  set  at  hot  end  while 
the  low-temperature  one  is  set  at  cold  end  in  accordance  with  the 
optimal  work  temperature  of  thermoelectric  materials.  The  simpli¬ 
fied  model  is  shown  in  Fig.  4.  The  performances  of  the  two-stage 
thermoelectric  unicouple  are  calculated  under  different  operating 
conditions. 

Figs.  5-7  show  that  the  junction  temperature,  output  power  and 
conversion  efficiency  of  the  two-stage  thermoelectric  unicouple 
rise  with  increasing  concentration  ratio.  When  the  concentration 
ratio  rises  to  173,  the  hot  end  temperature  of  the  low-temperature 
thermoelectric  unicouple  reaches  the  maximum  allowable  temper¬ 
ature  of  bismuth  telluride  of  about  225  °C,  and  the  concentration 
ratio  should  not  rise  any  more.  Meanwhile,  the  output  power 
and  conversion  efficiency  of  the  two-stage  unicouple  obtains  the 
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Fig.  3.  Dependence  of  the  output  power  on  temperature  difference. 


Fig.  6.  Dependence  of  the  output  power  of  the  two-stage  unicouple  on  concentra¬ 
tion  ratio. 
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Fig.  4.  Geometric  model  of  the  two-stage  thermoelectric  unicouple. 


Concentration  ratio 

Fig.  7.  Dependence  of  the  conversion  efficiencies  of  the  two-stage  unicouple  on 
concentration  ratio. 


Fig.  5.  Dependence  of  the  junction  temperature  of  the  two-stage  thermoelectric 
module  on  concentration  ratio. 

maximum  value  of  about  0.0353  W  and  13.31%,  respectively.  The 
contribution  of  the  low-temperature  thermoelectric  unicouple  is 
about  0.0279  W,  while  that  of  the  medium-temperature  thermo¬ 
electric  unicouple  is  about  0.0074  W.  It  is  found  that  the 


conversion  efficiency  of  the  two-stage  module  is  better  than  that 
of  the  single-stage  modules  because  the  thermal-electric  conver¬ 
sion  occurs  twice  when  the  heat  energy  passing  through  the 
two-stage  thermoelectric  module. 

3.3.  Three-stage  thermoelectric  module 

It  can  be  seen  from  the  previous  analysis  of  the  two-stage  mod¬ 
ule  that,  when  the  low-temperature  thermoelectric  module 
reaches  the  maximum  allowable  temperature  of  225  °C,  the  hot 
end  temperature  of  the  medium-temperature  thermoelectric  mod¬ 
ule  is  362  °C,  far  from  the  maximum  allowable  temperature  of 
skutterudite  of  about  525  °C.  Therefore,  it  is  feasible  to  add  an  extra 
medium-temperature  module  on  the  hot  side  of  the  two-stage 
thermoelectric  module  to  make  use  of  the  generating  ability  of 
skutterudite  between  362  °C  and  525  °C.  The  analysis  model  of 
the  three-stage  thermoelectric  module  based  on  one  low-temper¬ 
ature  module  and  two  medium-temperature  modules  is  estab¬ 
lished.  The  simplified  model  is  shown  in  Fig.  8.  Performance 
analysis  under  different  concentration  ratio  is  carried  out. 

Fig.  9  shows  that  when  the  concentration  ratio  rises  up  to  179, 
the  hot  end  temperature  of  the  low-temperature  thermoelectric 
unicouple  reaches  the  maximum  allowable  temperature  of 
bismuth  telluride  (225  °C),  while  the  medium-temperature 
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Fig.  8.  Geometric  model  of  the  three-stage  thermoelectric  unicouple. 


Fig.  11.  Dependence  of  the  conversion  efficiencies  of  the  three-stage  unicouple  on 
concentration  ratio. 


Fig.  9.  Dependence  of  the  junction  temperature  of  the  three-stage  thermoelectric 
module  on  concentration  ratio. 


Fig.  10.  Dependence  of  the  output  power  of  the  three-stage  unicouple  on 
concentration  ratio. 


thermoelectric  unicouples  is  operated  between  225  °C  and  513  °C 
approaching  to  the  maximum  allowable  temperature  of  skutteru¬ 
dite  of  about  525  °C. 


Table  2 

Performance  comparison. 


Unicouple 

Single-stage  (Bi2Te3) 

Two-stage 

Three-stage 

Concentration  ratio 

168 

173 

179 

Input  power  (W) 

0.2579 

0.2652 

0.2749 

Output  power  (W) 

0.0279 

0.0353 

0.0452 

Efficiency  (%) 

10.82 

13.31 

16.44 

Figs.  10  and  11  indicate  that  the  output  power  and  conversion 
efficiency  of  the  three-stage  unicouple  obtains  the  maximum  value 
of  0.0452  W  and  16.44%,  respectively,  at  the  ultimate  concentra¬ 
tion  ratio.  The  contribution  of  the  low-temperature  thermoelectric 
unicouple  is  about  0.0279  W,  while  that  of  the  two  medium-tem¬ 
perature  unicouples  are  0.0099  W  and  0.0074  W,  respectively.  It 
can  be  seen  that,  the  performance  of  the  three-stage  unicouple  is 
better  than  that  of  the  two-stage  unicouple  (13.31%).  The  extra 
medium-temperature  module  brings  an  efficiency  raise  of  3.13%, 
a  relative  change  of  23.5%.  Therefore,  the  design  method  of 
three-stage  thermoelectric  module  can  take  full  advantage  of  the 
characteristics  of  each  thermoelectric  material  and  effectively  im¬ 
prove  the  performance  of  power  generation. 


3.4.  Performance  comparison 

The  optimal  performance  of  the  single-stage,  two-stage  and 
three-stage  unicouples  are  obtained  from  above  analysis  and  the 
comparison  of  the  performances  is  summarized  in  Table  2. 

It  is  found  that  the  maximum  conversion  efficiency  of  the  ther¬ 
moelectric  unicouple  gradually  increase  with  increasing  stage 
number.  Therefore,  high  performance  of  thermoelectric  module 
can  be  achieved  within  the  large  temperature  difference  by  using 
a  multi-stage  thermoelectric  module.  Under  a  reasonable  condition 
of  the  efficiency  of  the  optical  focusing  system  rjopt  =  0.8  and  the 
absorptivity  of  the  heat  collector  r\a  =  0.8,  if  the  heat  losses  and 
the  contact  resistance  are  neglected,  the  utilization  of  solar  energy 
of  solar  thermoelectric  generator  can  reach  10.52%. 


4.  Conclusion 

A  3-D  finite  element  model  for  thermoelectric  module  based  on 
low-temperature  thermoelectric  material  bismuth  telluride  and 
medium-temperature  thermoelectric  material  filled-skutterudite 
is  established  and  performance  analysis  under  different  operating 
conditions  is  carried  out.  Result  shows  that,  reasonable  thermal 
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design  of  thermoelectric  module  can  take  full  advantage  of  the 
characteristics  of  each  thermoelectric  material  and  effectively  im¬ 
prove  the  performance  of  power  generation.  The  total  conversion 
efficiency  for  solar  energy  of  solar  thermoelectric  generator  based 
on  three-stage  thermoelectric  module  could  achieve  at  10.52%. 
With  the  continuous  emergence  of  new  thermoelectric  materials, 
solar  thermoelectric  generator  using  multi-stage  thermoelectric 
module  will  have  a  good  application  prospects. 
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